Renormalization group analysis of the spin-gap phase in the one-dimensional t-J model 
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We study the spin-gap phase in the one-dimensional t-J model, assuming that it is caused by the 
backward scattering process. Based on the renormalization group analysis and symmetry, we can 
determine the transition point between the Tomonaga-Luttinger liquid and the spin-gap phases, by 
the level crossing of the singlet and the triplet excitations. In contrast to the previous works, the 
obtained spin-gap region is unexpectedly large. We also check that the universality class of the 
transition belongs to the k = 1 SU(2) Wess-Zumino-Witten model. 
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The existence of a gap in the spin excitation has been 
considered to be a key to understand high-T c supercon- 
ductivity. This stimulated the study of one-dimensional 
(ID) electron systems some years ago. Recently, pos- 
sibilities of superconductivity in quasi ID systems are 
suggested 0j, and understanding of spin-gap phase in 
(quasi-) ID systems increases the importance. Now, we 
reconsider this problem from the ID t-J model which is 
the simplest but not fully understood. 

The Hamiltonian of the ID t-J model is written as 

ft = -t^ic^Ci+lc + 4+lv c i<r) 
i<j 

+J }XSj ■ Si+i - UfTii+i/4), (1) 



the spin gap generation mechanism is due to the attrac- 
tive backward scattering (scattering between electrons 
with the opposite momentum (kp,—kp) and spin) 
In this case, the universality class of the transition is the 
k = 1 SU(2) Wess-Zumino-Witten (WZW) model [§. 
On the basis of this assumption, we determine the tran- 
sition point with the singlet-triplet level crossing method 
P,ppC|] and we obtain the phase diagram (FIGjjj). Then 
we will verify the consistency of our method, considering 
the ratio of the logarithmic correction term. 

In general, the low-energy behavior of a ID electron 
system is described by the U(l) Gaussian model (charge 
part) and the SU(2) sine-Gordon model (spin part) 



;harge 
ill, 



H = H c + Ti s 



in the subspace without double occupancy. Generally, 
ID electron systems belong to the universality class of 
Tomonaga-Luttinger (TL) liquid which is character- 
ized by gapless charge and spin excitations and power-law 
decay of correlation functions. The phase diagram of the 
ID t-J model is obtained by Ogata et al, using exact 
diagonalization H . They found the enhancement of the 
superconducting correlation [K c > 1) and the phase sep- 
aration (K c — > oo) for large J/t region. They also found 
a phase of singlet bound electron pairs in the very low 
density region, but could get no evidence for a spin-gap 
phase by using a finite size scaling method at 1/3 fill- 
ing. Hellberg and Mele studied this phase by using a 
Jastrow-type variational wave function |^]. In their ap- 
proach, the variational parameter v is related with K c 
as K c = l/{2v + 1). They found that there exists a fi- 
nite region where the optimized parameter takes constant 
value v = —1/2 between TL phase and phase-separated 
state, and they interpreted the region as the spin-gap 
phase. Other variational wave function is proposed by 
Chen and Lee || . 

However, these authors did not discuss the detailed 
mechanism of the spin gap generation. One candidate of 
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Here a is a short-distance cutoff, g\ is the backward scat- 
tering amplitude, and for v = c,s 
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where n„ is the momentum density conjugate to <p v , 
[4> u (x), n„(a;')] = iS(x — x'), K v is the Gaussian coupling, 
and v c and v s are charge and spin velocities, respectively. 
The primary field of this model is exjr>i^/2(m u 4> v +n,j9 v ), 
where the dual field is defined as d x 9 u — 7rn„. In TL 
phase (gi > 0), the parameters K s and g\ will be renor- 
malized as K* = 1 and g\ = 0, reflecting the SU(2) 
symmetry. 

First, let us consider the case without renormalization, 
gi = 0. The finite size correction of the energy and the 
momentum of (p|) are described by the conformal field 
theory (CFT) Jl2||l3| with c — 1, where the central charge 
c characterizes the universality class of the model. For 
the t-J model, c = 1 as shown rigorously at J/t = 2 
p7[ and numerically M. The combined use of the CFT 
and the Bethe ansatz result gives a description of the ID 
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electron systems |14| , |15| , |16| , [L7| . The ground state energy 
of the system under periodic boundary conditions is given 
by 



E Q (L) = Le —rrj ~c, 

bL 



(4) 



where L is the system size. The excitation energy and 
momentum are related with exponents as 



2nv c 2nv s 
E — E = — —x c + —£-x Sl 

2tt 



(5) 



P - Po = -r{s c + s s ) + AkpD c + 2k F D s , (6) 

Li 

where kp — irN/2L with electron number N, and the 
scaling dimensions and the conformal spins are defined 
by x v = A+ + A~, s v = A+ — A~ , respectively, with the 
conformal weights, 




(7) 



The variables m„ and n v are related with electron quan- 
tum numbers as m c = 2D C + D s ,n c = AN c /2,m s — 
D s ,n s = AN S - AN c /2. Here AA C is the change of the 
total number of electrons, and AN S is the change of the 
number of down spins. D c (D s ) denotes the number of 
particles moved from the left charge (spin) Fermi point 
to the right one. N^ 1 (N^) is characterized by simple 
particle-hole excitations near right or left charge (spin) 
Fermi points. 

These quantum numbers are restricted by the selection 
rule under periodic boundary conditions [jlij 



AN C + AN, , 
D c = — n (mod 1), 



AN r 



(mod 1). 



(8a) 
(8b) 



In the case of twisted boundary conditions c 



t 

j+L,(T 



Cj a which is equivalent to the system where the flux $ 
penetrates the ring |TJ|, D c is modified as D c +&/2tt. For 
the ground state Eq, we choose periodic boundary condi- 
tions (<J> = 0) for N — Am + 2 electrons and antiperiodic 
boundary conditions ($ = 7r) for N = 4m electrons with 
an integer m. Changing the boundary conditions, the 
ground state becomes always singlet with zero momen- 
turn (P = 0) @§. 

In order to eliminate the contribution of the charge 
part, and extract the singlet and the triplet excita- 
tion in the spin part (x s — 1/2), we turn our at- 
tention on following states: (AN C , AN S , D c , D s ) = 
(0, ±1, 0, 0), (0, 0, Tl/2, ±1) under twisted boundary con- 
ditions ($ = 7r for N — 4m + 2, $ = for N = 4m). We 
can identify these excitation spectra by using (||) and 
(0), but the momentum P and the wave number p are 



not always identical. There is a relation P — p — &N/ L 
between them [^0| . 

Next, we consider the renormalization (g± 0). By the 
change of the cut off a — > e dl a, the coupling constant gi 
and are renormalizcd as M] 



dypjl) 

dl 
dyijl) 

dl 



-2/i 2 (0, 

-vo(l)yi(l), 



(9a) 
(9b) 



where yi(l) = gx/-KV s ,K s = 1 + y {l)/2. For the SU(2) 
symmetric case yo(l) = Vi(l), and yo(l) > 0, the scaling 
dimensions of the operators for singlet and triplet excita- 
tions \/2 cos v^^s (x ss ), and \[2 sin \/2</> s , exp(=Fiv / 2^ s ) 
(x s t) split logarithmically by the marginally irrelevant 
coupling as p2| 



yo 



1 3 

2 + Ay Q logL+1' 



_ l l 

^"^"I'wjlogL+l' 



(10a) 
(10b) 



where yo is the bare coupling, and we have set / = log L. 
This result is equivalent with that of the k = 1 SU(2) 
WZW model f|. Note that the ratio of the logarith- 
mic corrections are given by Clebsch-Gordan coefficients. 
When yo < 0, ya(l) is renormalized to yo(l) — * — oo, and 
there appears spin gap. At the critical point (ya = 0), 
there are no logarithmic corrections in the excitation 
gaps. The physical meaning of this point is that the 
backward scattering coupling changes from repulsive to 
attractive. And the SU(2) symmetry is enhanced at the 
critical point to the chiral SU(2)xSU(2) symmetry ||, 
since the spin degrees of freedom of the right and the left 
Fermi points become independent. Therefore, the critical 
point is obtained from the intersection of the singlet and 
the triplet excitation spectra d]|,[l(j ■ Using this method, 
we can determine the critical point with high precision 
|iof , since the remaining correction is only x s = 4 irrele- 
vant fields HQ. 

Here we analyze the numerical results for the t-J model 
(|l|) with the above explained method. We diagonalize 
L = 8-30 systems by the use of Lanczos and Householder 
method. An example of data [L = 16, n = N/L = 1/2) 
is shown in FIG.|[ Since the critical point is almost in- 
dependent of the system size as is shown in FIG]|, the 
phase diagram can be constructed without extrapolation. 
Our result is similar to the Hellberg and Mele's in the 
low density region, but the spin-gap phase spreads exten- 
sively toward the high density region. We are not able to 
answer whether the spin gap survives in the n — > 1 limit 
or not, because the numerical results become unstable 
in the high density region where the phase boundary is 
close to the phase-separated state. In TL phase, singlet 
and triplet superconducting correlations (SS, TS) have 
the same critical exponent 1/K C |1 |, while with a 
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spin gap, TS decays exponentially and SS is enhanced as 
l/K c , so that SS is dominant in the spin-gap region. 

In order to check the consistency of our argument, we 
calculate the ratios of the logarithmic corrections and 
scaling dimensions for the singlet and the triplet excita- 
tions from (|) and ([l0|). Here the spin wave velocity is 
given by p5 

E(L,N,S=1,P = 2*/L)-E (L,N) 

L^oo Z7T/ L 

which is extrapolated by the function v s (L) = v s (oo) + 
A/L 2 + B/L 4 . These corrections are explained by the 
irrelevant fields. The average of the renormalized scaling 
dimension (x ss + 3x st )/4, eliminating logarithmic correc- 
tions, and its finite size effect are shown in FIG.[| and 
FIG.[| respectively. The extrapolated data become 1/2 
with error less than 0.2 %. 

Finally, we discuss the reason why the previous studies 
have estimated the spin gap region very narrower than 
the real one. From the two-loop renormalization group 
equation of the k = 1 SU(2) WZW model ||]27|g| 



dyojl) 
dl 



the spin gap AE grows singularly as 



AE oc y/j - J c cxp(-Const./(J - J c )), 



(12) 



(13) 



where yo ot J c — J, therefore it is very difficult to find 
the critical point using conventional finite size scaling 
method. Note that ( |l3|) is the same asymptotic behavior 
as the spin gap of the negative U Hubbard model at 
half-filling, which can be obtained from the charge gap 
at positive U (2^], and the transformation between the 
charge and the spin degrees of freedoms |Q . 

In conclusion, we studied the spin-gap phase in the ID 
t-J model, considering the backward scattering effect in 
the TL liquid by the renormalization group analysis. Us- 
ing the twisted boundary conditions, we can extract the 
spin excitation spectra, and determine the critical point 
as in spin systems. The phase boundary is determined by 
the point where the backward scattering becomes repul- 
sive to attractive. The spin-gap phase obtained in this 
way is unexpectedly large, and the consistency of the ar- 
gument is also checked. This method can be applied to 
other models in ID electron systems, if the SU(2) sym- 
metry is assured. 
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FIG. 3. Size dependence of J c /t determined by the intersec- 
tions of the excitation spectra for L = 8, 12, 16, 20 systems at 
n = 1/2. These points are fitted by the form A+B / L 2 +C/ X 4 . 
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FIG. 1. Phase diagram of the ID t-J model (TL: TL 
phase, SG: spin-gap phase, PS: phase-separated state). In 
the spin-gap phase where the backward scattering is attrac- 
tive, the singlet excitation becomes lower than the triplet (see 
FIG.^|Q). The contour lines of K c are calculated by the data 
of L = 16 system |31|. 
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FIG. 4. Extrapolated value of (x S3 +3x 3 t)/4: and the scaling 
dimensions for the singlet (x ss ) and the triplet (x a t) excita- 
tions for L = 16 system at n = 1/2. 
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FIG. 2. Singlet and triplet excitation energies for L = 16 
system at n — 1/2. 
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FIG. 5. Size dependence of the averaged scaling dimension 
(x as + 3x 3t )/4 at n = 1/2. 
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